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We present here a new approach to determine an accurate variational wavefunction for general 
quantum antiferromagnets, completely defined by the requirement to reproduce the simple and well 
known spin- wave expansion. By this wavefunction, it is possible to obtain the correct behavior of the 
long distance correlation functions for the one dimensional S — 1/2 antiferromagnet, i.e. a system 
without long range order. The qualitative difference between the integer and half integer case is also 
easily understood with this variational approach. Finally we present numerical results for the 2D 
XY model, showing that the present wavefunction has an overlap > 0.99 to the exact ground state 
of the 5 = 1/2 model for finite system up to 6 x 6 clusters. 



PACS numbers: 75.10.Jm, 75.40. Cx, 75.40.Mg 

Since the discovery of High-Tc superconductivity an 
increasing attention has been devoted to the study of 
strongly correlated systems. In particular, the role of an- 
tiferromagnetic (AF) correlations in such electronic sys- 
tems was soon clear as it may lead to superconductivity 
at finite doping. |Q 

In the present paper we define a simple strategy to 
build a variational wavefunction for general quantum 
spin Hamiltonians, without restriction on dimensional- 
ity, spin, and existence or non existence of long range 
(AF) order. 

In order to simplify the discussion, we restrict the 
forthcoming analysis to the general anisotropic xxz 
model on a finite lattice with L sites and with periodic 
boundary conditions: 
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where S is usual notation for spin S operators and the 
couplings Jij , J z depend only on the distance between 
the sites i and j, in order to define a translation invari- 
ant Hamiltonian. It is also assumed that the couplings 
Js allow a stable ferromagnetic (FM) solution on the xy- 
easy plane, at the classical level. Within these notations, 
the spin isotropic Heisenberg model corresponds to near- 
est neighbor couplings with Jfj = Ji j for the AF case 
(after the usual transformation = (— iyS^ to change 
the sign of J in a bipartite lattice) and Jfj — — Jjy for 
the FM case. The isotropic condition is obtained by con- 
tinuity for | J z \ — > J. 

For this general Hamiltonian (|l|) standard spin-wave 
(SW) theory can be applied, where for convenience we set 
the order parameter along the y— axis and apply Holstein- 
Primakoff transformation rules at leading order in 
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Here the boson a] creates a spin fluctuation at the site i 
over the FM state \F > 

By introducing standard Fourier transformed variables 
a q , J(q) and J z (q), and after a little algebra, the linear 
SW Hamiltonian can be written in the compact form 
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where D q = 2J(0) - J(q) + J z (q) and r) q = -(J(q) + 
J z {q)). Each q component of the above Hamiltonian can 
be dealt independently and easily diagonalized by intro- 
ducing the Bogoliubov transformation a q = u q [3 q '+v q f3- q , 
with simple expressions for u q and v q in terms of D q and 
rjg. |H In a finite lattice system, a special care should 
be payed to the uniform q = mode corresponding to 
the conservation of the total spin component along the 
z— axis. As shown in the Bogoliubov transformation 
is singular for q — 0. Nevertheless, this mode can be dealt 
systematically, because it is exactly equivalent to a stan- 
dard projection of the SW ground state to the subspace 
of vanishing total spin projection S% ot — 0. Note that, 
in fact, the classical state \F > has not a definite S^ ot 
and after the above projection, indicated in the following 
by =o, the order parameter is uniformly distributed 
in the x — y plane, as is, of course, expected in a finite 
lattice. 

In view of the above discussion, the ground state wave- 
function has the following Gaussian form |0,|f : 
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The SW result is very simple because represents the so- 
lution of a quantum oscillator for S — > oo; however, in 
this limit the constraint, that the boson excitations on a 
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single site obey a\ati < 2S, is obviously violated in the 
wavefunction (||). 

We introduce here a simple wavefunction \tpr > which 
is for any S defined in the correct Hilbert space and only 
asymptotically for 5* ^ oo reduces to the known form 
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This is a generalization of the Manousakis wavefunction 
fiof obtained only for the isotropic AF model in 2D. The 
unknown function g q is then determined by requiring that 
the eigenvalue equation 
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is satisfied for S — > oo. Note that, due to the presence 
of the projector, the direction of the FM order param- 
eter in the x — y plane is irrelevant. It is thus con- 
venient to direct it in the x— axis direction before pro- 
jection. In the state \F >= FJt \Sf — S > each site 
has maximum spin along the x— axis, and in a basis 
diagonal with S z \a >= o~\o~ >, this state can be ex- 



panded as \Sf = S >= 2~ s 

m=-S 

therefore evident that the many-spin state Ps* = o I ^ > 
can be generally written as a sum of classical configu- 
rations |Cj >= | £7i, 02, • • • , &l >, with J^i °~i — with 
all positive, non vanishing coefficients. For instance, for 
S = 1/2, Ps~ ot =o\F > = sum all S z — spin configura- 
tions \C >. 

The exponential form in (^) (commuting with Ps^ ot ) 
represents a classical Jas- 

trow factor exp ^| ■ v(i — j)aia^ over the possible 

configurations |Cj > with a proper two-body potential 
v(r) = J] e~ zqR g q acting on the classical spins. The 

q^O 

positiveness of the wavefunction over all such configura- 
tions guarantees that (Q) is not orthogonal to the true 
ground state, for any finite S, and therefore it will nec- 
essarily collapse to it, provided (H) is verified. M 

We expect that the linear spin-wave approximation is 
quite accurate to determine the self-consistency in the 
eigenvalue equation (||), because this is determined only 
by the short range correlations, for which there is always 
some kind of order, needed to apply the basic approxima- 
tion (j^) . Provided the physics of the ground state wave- 
function is correctly described by (^), also long distance 
spin-spin correlations are expected to be correctly deter- 
mined. For instance, the spin- wave estimate is meaning- 
less for the order parameter m = S— < a^a >^ — oo 
for the ID S — 1/2 Heisenberg model, whereas with the 
present approach the long range order (LRO) of the wave- 
function Pgs —o\F > is readily destroyed by the long 
range potential v(r) ~ lnr, consistently determined (see 



later) in the large S limit. This approach allows therefore 
to define a spin-wave expansion even for models without 
a true LRO, which represents a remarkable extension of 
this powerful technique for quantum spin systems. || 

In order to determine the function g q , or equivalcntly 
to solve (|^) for S — > oo, we just notice that each q and 
—q couple of wavevectors in (|J) can be linearized by in- 
troducing complex auxiliary fields z q : 
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The above expression has to be applied to the vacuum 
of the spin waves a q \F >= 0. Since is, in this represen- 

one easily obtains that, for 



tation, S q — iy j (a q — 
large spin, 



-qj) 



9q 



9q 



t t 



By matching the two wavefunctions (^) and (|^), we de- 
termine g q — v q /(v q — u q ) for S — ► oo. 

The function g q , given by the above expression, turns 
singular only for q — > and behaving as ^ with 70 = 

— \fd for the nearest neighbor model. By expanding 
the exponential (^) in real space, we recover the wave- 
function from Ref. JtJ with logarithmic interaction be- 
tween the spins. As shown in this paper and related 
comments, this is precisely the condition to have Lut- 
tinger liquid behavior, or anomalous large distance expo- 
nents in ID. The spin-spin correlation functions, accord- 
ing to the Luttinger liquid analysis, decay as a power low: 
< SqS* >oc r~^ z and < SqS* >cx r ~^. Analogously to 
what was found in Ref. , we have obtained that fj, z and 
[i x depend only on the Jastrow coefficient 7: fi z = 
/j, x = — , in good agreement with the data shown in 
Fig. (jl]) . For the isotropic model the prefactor is not the 
exact one which is consistent with < So ■ S r >~ - corre- 
lations. However, if we determine g q by applying the SW 
expansion to the well known Haldane-Shastry model |Tl] ] 
( J(q) = \q 2 and J z (q) = —\{n-q) 2 mod(0,7r)), we find 
70 = —5, i.e. the exact value. The wavefunction \ipj- > 
is the exact ground state of the model for particular 
v(r) — 21nsin(^) to which our spin estimate is asymp- 
totically converging for L — * 00 ( g q = — + const.). 

Let us now consider the spin one model, again in ID. In 
our approach, the function g q is insensitive to the spin S, 
apart for a prefactor. However, the wavefunction changes 
dramatically, compared to the S — 1/2 case, due to the 
change of the Hilbert space, determined not only by two 
spins with opposite magnitude (S z = ±1), but also by 
vacancy sites with S z = 0, that are completely decoupled 
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from the non- vanishing spins. It is not difficult to realize 
that in the S = 1 case the model corresponds to an one 
dimensional Coulomb gas model with charge S z = ±1 
since the vacancy sites can be considered as an empty 
space between the charges. 
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FIG. 1. Log-log plot of the 5=1/2 spin-spin correlation 
functions c(L) for the spin operators in a plane (< Sq S~^^ 2 >) 
and in the direction orthogonal to it (< >) an d 

for two values of 7. The triangles refer to exact diagonal- 
ization data and squares to Monte Carlo estimates (error 
bars are much smaller then size of the points). The curve 
are obtained by fitting the data to the power law behavior 
log c(L) = /ilog L + a + b/L: 
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The Coulomb gas model with logarithmic interaction 
has been encountered several times in the literature and 
its phase diagram, probably complete, has been estab- 
lished. §]] But the mapping to this model is not ortho- 
dox, since with a little algebra the cos term in the Kane 
and Fisher Hamiltonian is replaced by a log(cos 2-S (<7<I>)) 
which differs from the usual model because there are infi- 
nite potential barriers between the valleys. However, we 
have found that the qualitative features of the cos model 
and the present logcos 25 one are the same. Note also 
that the mapping to a model with a local defect is possi- 
ble for integer S and not for half integer spin, where our 
wavefunction displays only one phase. 
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FIG. 2. Plot of the 5 = 1 spin-spin correlation 
function < 5g~5£/ 2 >> f° r Monte Carlo data, as a func- 
tion of 7. When 7 increase, the system is going from 
a disordered (plasma) phase, with a typical configuration 
+ H h - 0000 H h 000 - 00 + -0, to an ordered (dielec- 
tric) phase, with short range bound states with opposite spins 

— I 1 1 1 1 1 1 1 1 1 h The symbols 

+ ,0, — represent the values (5J = +1,0,-1). For the larger 
system size, the change of the phase is more pronounced, im- 
plying a true phase transition as expected (see text). 




FIG. 3. Log-log plots and L-log plots of the 5 = 1 in-plane 
spin-spin correlation function < 5q"5^ 2 > for 7 = 1.0 and 
7 = 1.8 (for the Monte Carlo data, with 7 = 1.0 error bars 
are smaller then the size of points). In the disordered phase 
(7 = 1.0), the spin-spin correlation function clearly decays 
as a power low, whereas the exponential decay assumption is 
much worse. 
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In the S = 1 case, for large 7, there is a dielectric phase 
(see Fig.||) that reminds the Haldane phase with an hid- 
den order parameter. Here the order is not hidden and 
of course the wavefunction is quantitatively meaningless. 
The remaining phase, where the charges are confined, is 
again characterized by a power law decay in the spin-spin 
correlations. This suggests that, for small 7, when our 
approximation is more reliable, a spin one magnet has 
power law correlations (Fig. |3|) and remains gapless. In- 
deed, the phase diagram of the S — 1 model [Q displays 
the power law behavior close to the J z ~ point, which is 
consistent with our findings as the estimated 7 decreases 
for smaller J z . In the Kane and Fisher model, correlation 
functions fall off with a power law in both regions. In the 
present model we cannot exclude the existence of a cor- 
rect finite correlation length (gap) for integer spins, but 
we have not been able to find a clear numerical evidence, 
as shown in Fig. (||) for 7 = 1.8. Nevertheless, we be- 
lieve that the emerging picture in the integer spin case is 
qualitatively sound and correctly describes the possible 
gapless phases, present also in the S = 1 models. 
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TABLE I. Difference between variational energy Ev and 
exact one E ex , and overlap square of the corresponding wave- 
functions for various clusters and number N of bosons. In 
this case the density of bosons was fixed by adding to the 
spin Hamiltonian (hi) a magnetic field term —HStot and de- 
termining the value of H by requiring: 2S(| — p) =< 5'°' > 
with p = j; fixed and consistently S — > 00. The direction of 
the order parameter has to be canted from the xy plane in 
order to determine a stable classical solution. 



Finally, we want to show a simple application of this 
variational wavefunction to a system of hard core bosons 
in the two dimensional case. Provided the system has 
LRO for p = N/L = 1/2 the spin- wave approxi- 
mation is expected to be accurate. Still, it was amazing 
for us to see so much accuracy (see table) in the ground 
state wavefunction, yielding practically an exact numer- 
ical solution of this correlated model. Notice that, at 
fixed number of bosons, the accuracy improves with size, 
showing that the low density limit || is fulfilled exactly, 
by this wavefunction. This is a non trivial fact because 
in this limit the Jastrow factor diverges logarithmically. 
Our numerical results shows that the spin-wave function 
(Q) is very accurate for any density, especially in the low 
density limit. We find therefore a clear evidence that the 
superfluid condensate exists, for any density, extending 
the rigorous proof in Ref. Jl2| for p = 1/2. 

To conclude we have presented here a very promising 
method to build unbiased variational wavefunctions for 
spin systems, which are quite reliable for gapless phases 
because we reproduce in this case all the known results 
for ID systems where the spin wave theory is worst. We 
believe that this wavefunction open new possibilities to 
understand strongly correlated systems. 
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